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ABSTRACT 

We explore the nucleosynthesis trends from two mechanisms during freeze-out ex- 
pansions in core-collapse supernovae. The first mechanism is related to the convection 
and instabilities within homogeneous stellar progenitor matter that is accreted through 
the supernova shock. The second mechanism is related to the impact of the supersonic 
wind termination shock (reverse shock) within the tumultuous inner regions of the 
ejecta above the proto-neutron star. Our results suggest that isotopes in the mass range 
12 < A < 122 which are produced during the freeze-out expansions may be classified 
in two families. The isotopes of the first family manifest a common mass fraction 
evolutionary profile, whose specific shape per isotope depends on the characteristic 
transition between two equilibrium states (equilibrium state transition) during each 
type of freeze-out expansion. The first family includes the majority of isotopes in this 
mass range. The second family is limited to magic nuclei and isotopes in their locality 
which do not sustain any transition, become nuclear flow hubs, and dominate the final 
composition. We use exponential and power-law adiabatic profiles to identify dynamic 
large-scale and small-scale equilibrium patterns among nuclear reactions. A reaction 
rate sensitivity study identifies those reactions which are crucial to the synthesis of 
radioactivities in the mass range of interest. In addition, we introduce non-monotonic 
parameterized profiles to probe the impact of the reverse shock and multi-dimensional 
explosion asymmetries on nucleosynthesis. Cases are shown where the non-monotonic 
profiles favor the production of radioactivities. Non-monotonic freeze-out profiles in- 
volve longer non-equilibrium nucleosynthesis intervals compared to the exponential 
and power-law profiles, resulting in mass fraction trends and yield distributions which 
may not be achieved by the monotonic freeze-out profiles. 

Subject headings: hydrodynamics — reverse shock — abundances, nucleosynthesis, 
radioactivities, supernovae: general — nuclear reactions 
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1. INTRODUCTION 



Despite the progress in core-collapse supernova theory during the past 20 years, the details 
of a fully self-consistent explosion and the nucleosynthesis are not yet fully understood. Three- 
dimensional supernova simulations with energy-dependent neutrino transport are necessary to clar- 
ify the growth of convective and Rayleigh-Taylor instabilities, the possibility of non-radial, non- 



axisymmetric instability modes, and the development of local fluid vortices (Janka et al. 2007) 
Further mixing processes near the proto-neutron star surface that could enhance the neutrino heat 



ing behind the shock include doubly diff"usive instabilities (Bruenn et al. 2004), neutrino-bubble 
instabilities ( [Socrates et al.||2005] ), and magnetic buoyancy instabilities ( [Wilson et al.||2005] ). The 
neutrino heating results in the outflow of baryonic matter (wind) from the surface of the proto- 
neutron star ( Duncan et al.[[1986[ ). This neutrino-driven wind interacts with the more slowly mov- 
ing, earlier supernova ejecta forming a wind termination shock, which is also termed as reverse 
shock ( Burrows et al.[[1995 ; Janka & Muller[[T995 ). The position of the reverse shock depends on 
the wind velocity and the thermodynamic conditions of the supernova ejecta which the wind col- 
hdes with ( [Arcones et al.[[2007[ ), while its boundary is strongly deformed to a non-spherical shape 
in multi-dimensional simulations ( [Arcones & Janka|2011| ). 



The interaction of the supernova shock with the stellar progenitor layers and the reverse shock 
with the mass layers above the proto-neutron star results in sudden entropy increases which trigger 
nucleosynthesis in the ashes of pre-coUapse nuclear burning. The expanding material then cools 
down until nuclear reactions freeze-out. The temperature and density trajectories of the expanding 
ejecta located initially outside the range of convective and mixing processes are monotonically 
decreasing. However, certain thermodynamic trajectories may be non-monotonic and can afl'ect 
the abundances' evolution in non-trivial ways. Possible factors to result in non-monotonic trajec- 



tories include the explosion energetics (Nakamura et al. 2001), explosion asymmetries in multi- 
dimensional simulations ( jNagataki et al.|199"7|[1998[ ), rotating progenitors ( [Fryer & Heger[[2000[ ), 
double explosion scenaria (Fryer et al. 2006[ ), and the eff"ect of the reverse shock (Arcones et al. 



2007, Arcones & Janka 2011). The yields from non-monotonic trajectories are usually diff"erent 
compared to compositions from monotonic trajectories. 

Certain studies have explored the impact of the supernova shock and the reverse shock on 
nucleosynthesis processes such as the a-rich and neutron-rich freeze-out ( [Woosley & Hoffman 
T992t [Witti et al.[[T994| ), the vjc-process ( [Frohhch et"aL][2006| [Pruet et al.[[2006l ), the rp-process 
( |Wanajo[[2006| ), and the r-process ( [Takahashi et al.[[T994l ). For the r-process specifically, the im- 
pact of the reverse shock has been studied under various wind schemes which include subsonic 



"breeze" solutions with constant pressure boundaries ( [Terasawa et al.[[2002[ ), supersonic winds 



with fixed asymptotic temperature ( [Wanajo et al.[[2002| ), and two-phase outflow models ( [Kuroda 
et al.[[2008] ), including reverse shock temperature ranges which could trigger the "cold" r-process 
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(|Wanajo|[2007|) or the classical r-process (|Kratz et aL]|2007[). However, integrated nucleosynthe 



sis based on those simulations shows that no heavy r-process elements can be produced ( Arcones 



& Montes 2011 1. In addition, reaction rate sensitivity studies have been performed to constrain 
uncertainties related to the p-process ( |Rapp et aL]|2006| ), vp-process ( jWanajo et al.||2011] ), freeze- 
out expansions from nuclear statistical equilibrium (henceforth NSE) or quasi-static equilibrium 
(henceforth QSE) in [The et al.| ( [T998] ); [Hoffman et aL| ( |20T0l ); [Magkotsios eTaLl ( |20T0l ), and hy- 
drostatic burning ( |Tur et al.]|2010[ ). Sensitivity studies add detail to the understanding of the mi- 
croscopic mechanisms which drive the related processes, because these mechanisms may not be 
easily identified for nominal values of the reaction rates or other related parameters. 

Yields of radioactivities as inferred from observations and presolar signatures (calcium and 
aluminum-rich inclusions, inter-planetary dust particles, and stellar grains) add complexity and 
raise new challenges for the explosion mechanism. Isotopic anomalies in presolar signatures indi- 
cate discrepancies between measured or inferred isotopic ratios and current supernova theory. The 
gamma-ray diffusive emission pattern of ^^Al indicates that massive stars are its most probable 
source ( [Leising & Diehl|[2009| [Limongi & Chieffi||2006bt |Tur et al.||2010[ ), given its correlation to 
^"Fe (Timmes et al. 1995[ Limongi & ChiefR 2006a). Yet, the synthesis of ^°Fe is largely dependent 
on the reaction rates producing it (Leising & Diehl 2009 1. "^'Ca found in SiC grains is likely to have 
common origin with ^^Al (ISahij pal et al.||l998 ), with a small fraction of '^^Ca possibly originating 
from supernova explosions (Nittler et al.|2008 1. It is argued that ^^Mn was uniformly distributed in 
the early solar nebula ( Yamashita et al.|2010 1. Although ^^Mn is produced in significant amount in 
core-collapse supernovae, it does not emit y-rays which would be easily detectable, but only X-rays 



which require significant accumulation in the interstellar medium for successful detection (Leising 
& Diehl|2009 ). '^Se also lacks a decay scheme by y-rays, although its detection in presolar grains 
is possible. Heavier isotopes such as ^^Zr ( |Lugaro et al.||2003 ), ^^Nb, and ^^Tc are believed to be 
produced primarily in asymptotic giant branch (AGB) stars, although we demonstrate cases where 
they may be produced in significant amounts during supernova explosions. 

In this study, we use monotonic adiabatic profiles to quantify nucleosynthesis trends caused 
by equilibrium patterns among nuclear reactions. In addition, we introduce non-monotonic param- 
eterized expansion profiles to simulate the impact of two mechanisms which may result in non- 
monotonic evolution for the local temperature and density evolution. The first mechanism involves 
multi-dimensional effects of the explosion following the passage of the supernova shock through 
the homogeneous layers of progenitor gas, and the second mechanism involves the impact of the 
abrupt deceleration by the reverse shock to the supersonic wind originating from the surface of the 
proto-neutron star. The mechanisms are independent of each other, because a contact discontinuity 
prevents mixing between the accumulated wind matter and the dense layer of shock-accelerated 



progenitor gas (Arcones et al. 2007). The non-monotonic profiles result in yield profiles which 
may not be achieved by the monotonic expansions. These results indicate the importance of multi- 
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dimensional supernova simulations on nucleosynthesis within the tumultuous stellar mass layers. 
We demonstrate cases where the non-monotonic profiles favor the production of radioactivities, 
and conclude that under certain circumstances non-monotonic profiles may increase the likelihood 
of detecting radioactivities in observed supernovae. 



Section[2]is a summary of our previous work ( [Magkotsios et al.|2010[ ). In Section|3]we present 



the formalism of the parameterized thermodynamic trajectories. Section|4]considers general trends 
of isotopes in the mass range 12 < A < 122 from our exponential and power-law trajectories. Sen- 
sitivities to reaction rate values related to the synthesis of radioactivities are discussed in Section|5] 
The impact of non-monotonic profiles on nucleosynthesis is discussed in Section [6} We conclude 
our discussion with a summary of our new results in Section |7] 

Our nomenclature is based on the following conventions. The temperature of the material 
is expressed by the pure number = T K). The electron fraction, or the total proton to 
nucleon ratio is Ye = Y^ZiYi = YjiZi/Aj X,. We define "nuclear flow" to mean the instantaneous 
rate of change of isotope i's molar abundance with time (dYi/dt) due to a given nuclear reaction 



(Iliadis 2007). For any single reaction linking isotope i with isotope j there is a forward flow, a 



reverse flow, and a relative net flow (/i;=(forward - reverse)/max(forward,reverse) that measures 
the equilibrium state of the reaction. 



2. SYNOPSIS OF PREVIOUS WORK 



The current work is a continuation of Magkotsios et al. (2010). We provide a synopsis of 



key notions and results presented in our previous work which are relevant to our new results. 



Magkotsios et al.| ( |2010[ ) used two parameterized expansion profiles within a parameter space of 
peak temperatures, peak densities, and initial electron fraction values, to discuss nucleosynthesis 
trends related to "^"^Ti and ^^Ni during freeze-out expansions from core-collapse supernovae. The 
initial composition was dominated by ^^Si, with neutrons or protons added to configure the value 
of the initial electron fraction Y^,. Setting the initial composition in this manner was not restrictive 
for the largest part of our parameter space. The peak conditions established a large-scale struc- 
ture composed by nuclear reactions in equilibrium and the mass fractions within this structure 
were rearranged to their equilibrium values early during the evolution. Large-scale equilibrium 
structures involve either nuclear statistical equilibrium (NSE), or global quasi-static equilibrium 
(global QSE). During NSE every reaction within a network is in equilibrium, while a small subset 
of them may be off equilibrium during QSE. For decreasing temperature and density values during 
the expansion the large-scale QSE cluster dissolves in multiple small-scale QSE clusters, where 
each cluster encompasses only a small number of nuclei with similar masses (local QSE). Further 
decrease in the temperature and density results in the cease of nuclear reactions (freeze-out stage). 
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Multiple types of freeze-out expansions were identified and discussed within the parameter 
space used, such as the normal freeze-out ( jWoosley et al.||1973| |Meyer|[T994| [Meyer et"aL]|1998 



Hix & Thielemann! 1999| ), the a-rich freeze-out ( [Woosley et al.||1973] ), the a-rich and proton-rich 



freeze-out (ap-hch), the incomplete silicon burning regime,, the photodisintegration regime the 
(p, y)-leakage regime for Yg > 0.5, the a-rich and neutron-rich freeze-out for Ye < 0.5 (orn-rich), 
and the depletion region of yields which separates the incomplete silicon burning and normal 
freeze-out regions from the or-rich freeze-out (chasm). 

Each type of freeze-out is related to distinct transitions between two very different equilibrium 
states. In this work we refer them as "equilibrium state transitions" and we abbreviate them as EST. 
These transitions between equilibrium states involve multiple reactions within the QSE cluster 
which break equilibrium, resulting in changes to the QSE cluster's size and shape. The scale of the 
transitions during freeze-out expansions may be global, such as the division of a large QSE cluster 
encompassing almost every isotope in the network into two large QSE clusters localized within 
the silicon and iron groups respectively. Alternatively, the scale of the transitions may be local, 
involving only a small set of nuclear reactions which form small-scale clusters such as equilibrium 
chains of (p, y) or (n, y) reactions interacting with a-captures or (p, n) reactions, respectively. ESTs 
are entropy driven, where the temperature sets an approximate threshold for a transition, while the 
density at the threshold temperature determines whether the transition takes place or not. Electron 
fraction variations, the expansion timescale, and key reaction rates control the local equilibrium 
patterns which shape the locus of each region. 



3. PARAMETERIZED PROFILES 

The nucleosynthesis calculations implement mature reaction network solvers ( Timmes|1999 



Fryxell et al.||200 0) and utilize the 489 and 3304 isotope networks described in |Magkotsios et al. 



( |2010 j). The 489 isotope network has been expanded to 553 isotopes to include a sufficient amount 
of isotopes near the magic number 50, and is listed in Table [T] The 553 isotope network includes 
two separate isomers for ^^Al, ^^Alg for the ground state and ^^Al^ for the first excited metastable 
state. We use three parameterized profiles to model the freeze-out expansions. All profiles assume 
that a passing supernova shock wave heats material to a peak temperature Tq and compresses the 
material to a peak density po- This material then expands and cools down (freezes out) until the 
temperature and density are reduced to the extent that nuclear reactions cease. 

The first two parameterized profiles involve monotonic temperature and density decrease, 
following a constant T^/p evolution which implies constant radiation entropy in suitable limits. 
The first profile is the exponential expansion ( Hoyle et al.|1964[ [Fowler & Hoyle|1964 ) 



Tit) = ToCxp{-t/3T) p(0 =Poexp(-f/r) (1) 
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with a static free-fall timescale for the expanding ejecta 



(2) 



The exponential profile has been used extensively in the past to explore yield trends and their 
sensitivity to reaction rates or electron fraction values ( [Woosley et al.||1973| [Woosley & Hoffman 
T992t [The et al.|[T998| [Hottman et aL][20T0t [Magkotsios eTaTHlOTOl). The second profile is a 



power-law based on homologous expansion introduced by [Magkotsios et aL] ( |2008[[2010[ ) 

Tq , Po 



T(t) = 



ct^-\ 



{ct+\f 



(3) 



where the coefficient c = 2 s ^ is chosen to mimic trajectories taken from core-collapse simulations. 

Figure [T] compares the general properties of the exponential and power-law profiles. For a 
given initial condition, the power-law evolution is always slower than the exponential one. More- 
over, the power-law evolution becomes slower for increasing peak temperature and density values. 
The differences in these two profiles affect the final yields as material traverses different burning 
regimes on different timescales. The figure also depicts the NSE, global QSE, local QSE, and final 
freeze-out burning regimes. The exponential and power-law trajectories are chosen so that they 
generally bound the temperature and density trajectories of hydrodynamic particles from spheri- 
cally symmetric and two-dimensional explosion models. 

We introduce non-monotonic profiles to simulate possible thermodynamic conditions in the 
ejecta from different mass layers within the star. Non-monotonic temperature and density evolu- 
tion may arise due to convection and instabilities following the heating of homogeneous progenitor 
mass by the supernova shock, or by the termination shock to the supersonic wind within the in- 
ner mass layers above the proto-neutron star. Our profiles involve three stages to model such 
possible effects (Figure [1]). During the first stage, the supernova shock rises the temperature and 
density values of the material traversed, and the material is allowed to expand while it cools and 
rarefies. During the second stage we introduce a simplified approach to model the effect of multi- 
dimensional asymmetries to the explosion or the reverse shock within a parameterized expansion 
profile. This stage involves a contraction phase which rises the temperature and density linearly in 
time to a local maximum. The third stage involves an exponential freeze-out, because the material 
is assumed to be part of the ejecta and eventually escape from the star. Contrary to the exponential 
and power-law cases, no explicit assumption is made about holding Ip constant. 

Each temperature and density trajectory has a peak value followed by a local minimum and 
then a local maximum. The local minimum and maximum values for the case of the reverse 
shock depend on the wind velocity and the deformed boundary of the shock near the supernova 
ejecta which the wind collides with ( Arcones et al.|20()7 ; Arcones & Janka||2011| ). For the case of 
our parameterized expansion profiles, we choose the local extremum values for temperature and 
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density randomly from a uniform distribution. Each value of a local extremum point is considered 
to be independent and identically distributed. This is a reasonable approach, because we aim 
to study the key nucleosynthesis trends within the tumultuous inner layers of the ejecta, and we 
attempt to simplify the hydrodynamic evolution for this purpose. 



The differential equations for the non-monotonic temperature and density trajectories are 



dT 

dt 
dp 
dt 



-Toexp(-k^t) 
-po exp(-kpt) 



k^b^t + - — — + - — ^—r 
\ +a^t (1 +a^0 



kpbpt + ^ + -- — ^— T - bp 



and their solutions are 



r(0 = roexp(-A:,0 
pit) = Po exp(-kpt) 



+b^t 

\ +aj 



+ bj 



1 + apt 



(4) 
(5) 

(6) 
(7) 



where the parameters a,, bj and kj may be chosen to control the local extremum points. For 
instance, the density minimum is given approximately by p„„>, = po(2^Ja^ - bp) lap at time 
tmin = ( - ^Jbp)|{ap ^/Fp), while the maximum is given by p,r,ax = Po/(e kp) at time t^^x = 1 /kp. 
The red curves of Figure [T] show the profile's general trends compared to the exponential and 
power-law expansions. The ascending trajectory is focused on QSE and local equilibrium stages. 
We ensure that the local maximum for the temperature does not exceed the NSE threshold, other- 
wise the preceding part of the trajectory would not impact the evolution. The local minimum for 
the temperature ranges from the freeze-out temperature Tg = 0.01 until the peak temperature. The 
subsequent local maximum for the temperature ranges from the value of the local minimum until 
the value = 4. Thus, a non-monotonic evolution is guaranteed without the re-establishment of 
NSE. The range of the extremum point values for the density profile is less constrained and even 
monotonic profiles may arise. The minimum spans the range 10^ g cm"^ < p„j,„ < po and the range 
for the maximum is 10^ g cm"^ < p^^x < 10^ g cm~^. 



4. YIELDS FROM THE EXPONENTIAL AND POWER-LAW PROFILES 

The exponential and the power-law expansion profiles used in this work set the basis for 
quantifying the details of nucleosynthesis mechanisms. Their low computational cost allows the 
monitoring of large-scale and small-scale equilibrium patterns among nuclear reactions in time, 
which is a powerful tool for identifying microscopic components that affect the composition of the 
ejecta. 
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Monotonic profiles are probes for multiple burning regimes. For instance, during subsonic 
outflows (neutrino-driven "breezes") from the proto-neutron star, the flow merges smoothly with 
the denser shell of ejecta behind the outgoing supernova shock resulting in monotonic evolution 
for both the temperature and density of the ejecta ( jOtsuki et al.||2000[ |Terasawa et al.|[2002| |Ar- 



cones et al.|2007 1. For cases where the temperature increase imparted by the reverse shock during 
supersonic winds is above the NSE threshold T<) ~ 5, the nuclear abundances acquire NSE values 
and the previous thermodynamic history of the ejecta is irrelevant to the nucleosynthesis evolution 
following the temperature jump. Temperature increases above the NSE threshold may also occur 
when the homogeneous progenitor matter is traversed by the supernova shock. Monotonic profiles 
are suitable probes for cases where mixing processes below the NSE threshold are either absent or 
negligible. 

Figure |2] shows the final mass fractions of isotopes along the or-chain for initially symmetric 
matter {Ye = 0.5). The temperature-density planes include the full range of peak conditions within 
our parameter space. With the exception of ^^Ni (not shown), the topological structure of all planes 
is similar, and is marked by distinct regions. These regions are labeled in the ^^Si contour plot for 
the symmetric case, each region corresponding to a type of freeze-out expansion. Further types 
of freeze-out expansions are manifested for initial electron fraction Yg 4^ 0.5 (see analysis below). 
The aggregate range of isotopes produced by all identified freeze-out types is in the mass range 
12 < A < 122, including in addition the free neutrons, protons, and or-particles. The structure of 
the temperature-density plane is very similar among the majority of the isotopes in this mass range. 
We classify into a family the isotopes within the mass range 12 < A < 122 whose temperature- 
density plane features a region-divided structure (henceforth the "first family" of isotopes) and 
explore the common features that these isotopes share. The or-chain isotopes shown in Figure [2] 
belong to the first family. 

The mass fraction profiles per region for the isotopes in Figure [2] and the first family overall 
are similar, indicating that within a region the isotopes of this family are produced by the same 
mechanism. The mass fraction profile similarities arise from the initial formation of a large-scale 
QSE cluster during the freeze-out processes. Within the cluster all nuclei are interconnected with 
reactions in equilibrium, and mass fraction values are determined by the temperature, density and 



electron fraction variations based on minimization principles of the Helmholtz free energy (Seiten 



zahl et al. 2008 ). Subsequent ESTs alter the shape of the QSE cluster and eventually the cluster 
dissolves. The precise locus of the regions in the temperature-density plane for an isotope of the 
first family depends on the local equilibrium patterns near the isotope while the large QSE cluster 
dissolves. 



Figure [3] illustrates the upward shifting in mass of the QSE cluster ( [Meyer et al.||1998 ). The 
QSE cluster remnant condenses around the magic number 28, and nuclei in this small group tend 



-9- 



to dominate the final composition. The isotopes of the first family which are gradually left outside 
the QSE cluster form chains of (p, y) reactions in equilibrium along the isotone lines. The first 
EST related to these isotopes' exit from the QSE cluster is signaled at the microscopic level by 
the equilibrium break of the a-capture reactions linking the (p,7) equilibrium chains. During the 
ap-rich freeze-out, isotopes of the first family sustain a second EST when certain (p, y) reactions 
in the isotone chain break equilibrium. These small-scale equilibrium patterns are responsible for 
producing eventually the isotopes of the first family from ^^C to the iron peak. On the contrary, 
the formation of the chasm for each isotope of the first family results from the dissolution of 
the large-scale QSE cluster to two smaller ones. The first cluster encompasses the silicon group 
elements and the second cluster encompasses the iron group elements. The cluster breakage results 
in massive flow transfer from the silicon and most of the iron group isotopes toward a small group 
of nuclides near the magic number 28. The flow transfer proceeds until all mass fractions are 
depleted, excluding the mass fractions of nuclei around the magic number 28. These nuclei are 
produced in large amounts and dominate the final composition. 

The types of freeze-out discussed so far (normal, a-rich, ap-rich, and the chasm) tend to favor 
the production of nuclei with proton and neutron numbers in the locality of the magic number 28. 
Figure |4] shows a sample of such nuclei, and Table |2] provides the complete list. These isotopes 
tend to dominate the final composition for most initial electron fraction values. The final mass 
fractions in Figure |4] demonstrate homogeneous structures within the temperature-density plane, 
implying that these isotopes do not sustain any EST during the evolution. The restriction of the 
remnant QSE cluster and the accumulation of nuclear flow among these isotopes are responsible 
for the absence of ESTs. The accumulation of flow stems from the fact that nuclei with proton or 
neutron numbers near the magic number 28 tend to maximize their binding energy per nucleon. 
As a result, such nuclei are relatively more bound compared to nuclei with nucleon numbers far 
from the magic number values, and their production within a network of reactions is favored. We 
classify isotopes which do not sustain any EST during freeze-out expansions and tend to dominate 
the final composition into a "second family" of isotopes. We have demonstrated that nuclei whose 
neutron or proton number is near the magic number 28 belong to the second family. Below, we 
show that nuclei with neutron numbers near the magic numbers 50 and 82 also belong to the second 
family. 

Figures |5] and [6] show the temperature-density planes of select radioactivities up to mass 
A = 97 which have non-negligible yields for the corresponding initial Ye values. The regions 
of the orn-rich freeze-out and (p, y) -leakage regime are labeled. These two types of freeze-out 
expansions are not manifested for initially symmetric matter. The temperature-density planes 
depict a region-divided structure, indicating that these radioactivities belong to the first family of 
isotopes. The decay timescale of ^^Al^ is approximately 2 s. Consequently, yield values are seen 
only for the exponential profile where the expansion timescale is less than a second, while for the 
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power-law it decays prior to complete freeze-out. On the contrary, ^^AL is mostly produced during 
the power-law expansion for > 0.5. '^'Ca is produced during the a-rich and ap-rich freeze-outs 
for the full range of our initial electron fraction values, and also during the (p, 7)-leakage regime 
(Figure [6]) for the exponential profile and Yg > 0.5. The (or, y) and (a,p) channels control its 
production for Y^ = 0.48, while the (p,7) and the weak reactions impact its synthesis for Y^ > 0.5. 
In addition, the (a, p) channels shape the locus of the borderline between the (p, 7)-leakage and 
Si-rich regimes in the contour plot for Y^ > 0.5. ^^Mn is produced mostly during the normal 
freeze-out for Y^ < 0.5, although it has significant yields from the a-rich and ap-rich freeze-outs 
for the power-law expansion. It is relatively insensitive to reaction rates for Ye < 0.5, while the 
(p, y) and weak reactions have significant impact for Ye > 0.5. The (a, p) and (p, y) channels affect 
the borderline between the (p, y)-leakage and Si-rich regimes. 

Radioactivities heavier than mass A = 60 are produced primarily in neutron-rich environments 
for the exponential and power- law profiles. Specifically, they may either be produced during an 
an-rich freeze-out ( |Woosley & Hoffman|1992 ) or by a process that combines features between the 



a-rich and an-hch freeze-outs (Figure [5]). The an -rich freeze-out occurs for relatively low initial 
electron fraction values and its locus is constrained to regions of low peak densities in the contour 
plots. The combination of high peak temperatures and low peak densities allows the establishment 
of a photodisintegration regime early in the evolution. The balance between the p(e",ye)n and 
n(e^, v^)p reactions maintains the electron fraction values below 0.5, which favor an overproduction 
of neutrons against protons. Such values for Ye allow the major nuclear flows to bypass the doubly 



magic nucleus ^^Ni and heavier elements are produced (Hartmann et al. 1985 Woosley & Hofi'man 



1992, Magkotsios et al. 2010). The QSE cluster shifts upward in mass, but it is not localized 
solely around the magic number 28. Instead, neutron capture reactions shift the cluster to heavier 
masses and pile up nuclear flow in the locality of nuclei with neutron magic numbers 50 and 82. 
The concentration of nuclear flow around these nuclei maintains the equilibrium structure in their 
locality and prevents them from sustaining ESTs. Once again, the flow concentration near nuclei 
with magic numbers 50 and 82 is a nuclear structure effect. These nuclei maximize locally the 
nuclear binding energy per nucleon and are relatively more bound compared to other nuclei with 
nucleon numbers away from the magic number series. Consequently, isotopes such as ^^Kr, *^'Rb, 
^^Sr, and ^^^Zr belong to the second family of isotopes and dominate the final composition along 
with free a-particles and neutrons. During the evolution, the excess of free neutrons guarantees 
a large-scale equilibrium structure maintained primarily by chains of (n, y) and (p, n) reactions in 
equilibrium. Isotopes of the first family with mass A > 60 including the radioactivities ^"Fe and 
^^Se sustain an EST when (p, n) reactions in their locality break equilibrium. 

The production of radioactivities such as ^^Nb, ^^Zr, and ^^Tc by the exponential and power- 
law expansions is favored only within the narrow transition region between the or-rich and orn-rich 
freeze-outs. Figure |7] shows the mass fraction evolution of dominant elements and radioactivities 
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within this region. The free neutrons are depleted below the NSE threshold and do not allow the 
neutron capture reactions to shift the QSE cluster until the magic number 82. The flows are blocked 
around the magic number 50. This effect is illustrated in Figure [sj The yields of **^Kr, ^^Rb, and 
^^Sr still dominate the final composition, but are slightly enhanced compared to the region of the 
an-rich freeze-out. The constraint of nuclear flow near ^^Kr maximizes the yields of ^°Fe and ^'^Se, 
and allows the production of ^^Nb, ^^Zr, and ^^Tc (Figure [s]). 



5. REACTION RATE SENSITIVITY STUDY FOR SELECT RADIOACTIVITIES 



We perform a sensitivity study on reaction rates related to the synthesis of radioactivities for 
the exponential and power-law expansions, to identify the reactions which are primarily responsi- 
ble for their production in core-collapse supemovae. These critical reactions determine whether an 
EST takes place or not. Reaction channels and individual rates are either multiplied by a factor or 
removed from the network. Strong reaction rates are multiplied by factors of 100 or 0.01, while the 
corresponding factors for the weak reactions are 1000 or 0.001. These factors are adequate to fa- 
cilitate the identification of trends in the yields, although they exceed experimental uncertainties in 
most cases. Sensitivity studies where reactions were varied within experimental uncertainty ranges 
have been performed by Hoffman et al. (2010) and Tur et al. ( 2010[ ). Our calculations include rates 
for weak interactions ( jFuUer et al.||1980[ |1982b|a[ |Oda et al.]|1994^ |Langanke & Martinez -Pinedo 



2001), the theoretical rates of Rauscher & Thielemann (2000), and select experimental rates for 
capture and photodisintegration reactions. 

The structure of the temperature-density plane for isotopes of the first family is affected by 
certain key reactions per isotope, in combination with the expansion timescale. Specific reactions 
such as the 3a, aan, p(e , Ve)n, n(e^, v^)p, and combinations of a large number of weak reactions 
impact all mass fractions simultaneously either by transferring nuclear flow to the QSE cluster, or 
by contributing to electron fraction variations ( [Fuller & Meyer] 1995^ [McLaughhn & Fuller||1995 



McLaughlin et al.|1996t|Aprahamian et al.|2005HSurman & McLaughlin|2005| [Liebendorfer et al 



2008[ ). However, the local equilibrium patterns are controlled by reactions in the locality of each 
isotope. Table [3] lists the reactions that impact the synthesis of radioactivities in the mass range 
12 < A < 122 produced during freeze-out expansions. In Table [3] the contribution of a reaction 
is focused on specific parts of the isotope's mass fraction curve. Below we use the term "arc" 
frequently, so it is convenient to provide a visualization of this structure. Figure |9] shows a sample 
mass fraction evolution of an isotope for decreasing temperature. Two local minimum points of 
the mass fraction curve are identified. These points separate the curve in three parts (arcs). The 
first part (black arc) is related to the mass fraction evolution when the isotope participates in a 
large-scale QSE cluster (global QSE). The second and third parts (red and blue arcs respectively) 
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are related to mass fraction trends when the isotope either participates in small-scale equilibrium 
clusters (local QSE) or does not belong to any cluster at all (non-equilibrium nucleosynthesis). 
Since the third arc is not a full arc, it is also mentioned as "ascending track". Note that the mass 
fraction curve may be limited to two arcs, i.e. the first arc and the ascending track until freeze-out 
(for instance, see Figure [7]). 

The ^^Al yield has an average value X(^^Al) ~ 10"^ within the ar-rich and ap-rich freeze-out 
regimes for the exponential profile only, where the metastable state ^^Alm has significant yield 
(Figures [5] and |6]). ^^Alg(p,y)^'Si and ^^Alg(ff, p)^^Si configure the chasm features and the char- 
acteristic arcs in the mass fraction profiles during the or-rich and ap-rich freeze-outs. They also 
control the flow between the A'^ = 13 and A'^ = 15 isotones. The equilibrium break of ^^Alg(Qr, p)^^Si 
marks the appearance and controls the depth of the chasm in the contour plot and the first dip in the 
^^Al mass fraction profile. The equilibrium break of ^^Alg(p,7)^'Si during the ap-rich freeze-out 
transfers flow from ^^Al to heavier isotopes along the = 13 isotone and configures the shape 
and dip of the second arc in the ^^Al mass fraction. Additional reactions that contribute similarly 
are ^^Mgip, yf^Alg, ^^AKp, yf^Si, ^^Si(p, yf^F, and ^'^Fip, yf^S. The '-^Mg(a, yf^Si reaction dis- 
tributes flow within the silicon group and shapes the ascending track to the ^^Al mass fraction past 
the arcs during the ap-rich freeze-out. The yield of ^^Al is largely dependent on the collective 
flow transfer by weak reactions toward symmetric nuclei. Specific weak reactions with the largest 
contribution for ^^Al are listed in Table lU 

"^'Ca is an example of composite contribution from (p, 7) reactions along its own and neigh- 
boring isotone lines which are connected by (a, 7) reactions. Table[3]lists the related reactions for 
Ye > 0.5, and reactions that transfer flow for neutron-rich compositions. For proton-rich environ- 
ments ^^Ca is also dependent on the collective flow transfer by weak reactions. 

Radioactivities with mass A > 60 depend mostly on neutron captures and weak reactions. The 
yield of ^''Mn depends only on the collective flow transfer by weak reactions. The weak reactions 
with the largest contribution are listed in Table [3] *^^Kr(a, n)^^Sr is the main flow distributor within 
the small cluster in the locality of the neutron magic number N = 50 and it affects the mass 
fractions of ^"Fe, ^^Se, ^^Nb, ^^Zr, and ^^Tc. ^°Fe(p, n)^''Co shapes the second arc of the ^°Fe mass 
fraction, and its strength determines the degree of the yield's depletion. ^*'Fe(a, n)^^Ni contributes 
to the formation of the arc, while ^^Fe(n, 7)^°Fe and ^^Fe(a, n)^^Ni regulate the arc's amplitude and 
slope respectively. Most notably, ^°Fe(n, 7)^^Fe does not impact the ^°Fe mass fraction and yield for 
the freeze-out expansions. The mass fraction profile of ^^Se is marked by a sharp ascending track at 
complete freeze-out which increases the yield by an order of magnitude. The abrupt flow transfer 
stems mostly from the collective contribution of the weak reactions, with major contribution from 
^^As(, e~ Ve)^^Se. Additional reactions that contribute to this ascending track prior to freeze-out 
are '^Se(n,7)^^Se and ^^Se(n, 7)''^Se. Its mass fraction arc is formed by ^^As(p, n)^^Se. Similarly, 
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the mass fraction arc for ^^Nb is formed by ^^Zr(p, n)^^Nb, and the arc for ^^Tc by ^^Mo(p, n)^^Tc, 
^'Tc(n, 7)^^Tc, and ^^Mo(p,n)^^Tc. ^^Zr is in the locality of the neutron magic number N = 50, 
and its mass fraction profile monotonically increases up to X^^Zr) ~ 10"'^ in the transition region 
between the or-rich and an-rich freeze-outs. The rate strength of ^^Zr(n, yf^Tx shapes the yield 
value by transferring flow to ^^Zr from the isotopes with neutron number A'^ = 50. 



6. YIELDS FROM NON-MONOTONIC PROFILES 

The exponential and power-law profiles discussed so far involve the initial formation of a 
large-scale QSE cluster and subsequent states of small-scale clusters. Non-equilibrium nucleosyn- 
thesis appears only during the freeze-out stage at the end of the evolution. The non-monotonic 
expansion profiles (Section [3]) may have non-equilibrium intervals followed by the formation of 
a large-scale QSE cluster, resulting in final compositions which cannot be achieved by mono- 
tonic profiles. Table |4] lists the diff'erences between the monotonic and non-monotonic profiles. 
The keywords related to the evolution of the QSE cluster are "hierarchical" and "periodic" for 
the monotonic and non-monotonic profiles respectively. The term hierarchical denotes the grad- 
ual dissolution of the single large-scale QSE cluster to multiple small-scale clusters, while the 
term periodic describes the sequence of transitions from NSE and global QSE, to local QSE and 
non-equilibrium phase, to global QSE again and then local QSE and non-equilibrium phase until 
freeze-out. 

The existence of ESTs during non-monotonic profiles depends on the combination of the 
extremum point values for the temperature and density. Arcones et al. ( |2007 ) report ranges of 



0.4 < < 2 and 10^ < p < 10"^ g cm ^ for the temperature and density behind the reverse 
shock, while Wanajo et al. ( |2011| ) report a temperature range of 1.5 < < 3. The flow transfer 



by reactions out of equilibrium among scattered small-scale equilibrium clusters impacts the mass 
fraction evolution dramatically. The flow patterns are very sensitive to variations in the values of 
the temperature and density extremum points. This sensitivity diversifies the production mecha- 
nisms significantly. Table [5] lists the combinations of temperature and density extremum points 
within our data set which tend to maximize the yields of radioactivities. The peak temperature for 
the non-monotonic expansions is chosen to be = 9, and the peak density and initial electron 
fraction values range between 10^ < p < 10*^ g cm"^ and 0.48 < 7^ < 0.52 respectively. These 
peak temperature and density values are large enough to establish NSE early in the evolution, so 
that the initial composition dependence is removed. Below we present the details of certain pro- 
files which produce simultaneously most of the radioactivities in the mass range 12 < A < 122. 
These radioactive isotopes could be transferred to upper (cooler) mass layers by mixing processes 
during the explosion, where the eff'ective lifetime of the radioactivities is longer ( |Tur et al.|2010| ). 
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and possibly released to the interstellar medium. 



Figure 10 shows the thermodynamic and mass fraction evolutions of two non-monotonic pro- 
files for initially neutron-rich composition (Ye = 0.48). The extremum values for the profile of the 
top row are Tf" = 0.02, T:;'"' = 2.5, p„,in = 2 x 10^ g cm-\ and p,nax = 3 x 10^ g cm^l This profile 
approximately maximizes the ^"Fe yield within our data set. During the initial temperature decre- 
ment (until t = 10 s) the density is relatively fixed to its peak value. The conditions at the Tg = 5 
NSE threshold are similar to an cr-rich freeze-out in neutron-rich matter, where the 3a forward 
rate dominates its inverse photodisintegration and the protons are rapidly consumed. However, the 
temperature evolution is fast enough that it prevents significant flow transfer beyond '^O, and the 
QSE cluster dissolves without shifting upward in mass. The temperature nearly reaches freeze-out 
levels after t = 10""^ s, and ^^C dominates the composition, with significant mass fraction values for 
free neutrons, a-particles, and '^O. The density decrease from t = 10'"^ s until t = 10~^ s does not 
impact the mass fractions due to low temperature values. The subsequent temperature and density 
increase result in carbon burning primarily by the '^C('^C,p)^^Na and '^C(^^C, a)^*'Ne reactions. 
A second large-scale QSE cluster is formed and the low electron fraction values in combination 
with the free neutron abundance guarantee flow transfer near ^^Kr through (n, y) and (p, n) reac- 
tions. The temperature maximum is slightly lower than the typical silicon burning threshold ~ 3 
( Iliadis|2007 1 at ? = 1 s, and the final composition is dominated mostly by carbon burning products 



and elements in the locality of the A'^ = 50 neutron magic number, while elements near the N = 2S 
magic number are severely underproduced. The final composition includes significant yields for 
the radioactivities listed in Table [sj excluding '^'^Ti and ^'Tc. 



The second row of Figure|10|corresponds to a profile with extremum values 1'™" = 1 .4, Tg""^ = 
1.94, pm„, = 3.3 X 10^ g cm"^, and p,„ax = 3.4 x 10^ g cm"^. This profile approximately maximizes 
the '*'*Ti yield within our data set. The flow transfer by the 3a rate to the QSE cluster occurs 



over a longer interval compared to the profile of the top row in Figure 10 The QSE cluster shifts 



upward in mass until the N = 50 neutron magic number, producing Kr and the radioactivities 
^^Nb, ^^Zr and ^^Tc. Up to the point where the ^^Kr abundance is maximized (near t = 10"^ 
s), the process resembles the evolution of exponential and power-law profiles in the transition 
regime between the a-rich and an-rich freeze-outs (see Section |4] and Figure |7]). The first a- 
capture reactions to break equilibrium appear in the mass region of ^^C and '^O. The increasing 
number of a-capture reactions which break equilibrium within this mass region and the silicon 
group shift the low mass border of the QSE cluster toward heavier nuclei. The a-capture reactions 
with the largest net flows are blocked within the silicon group, because the Coulomb repulsion for 
the specific thermodynamic conditions prevents the a-capture reactions involving heavier nuclei 
from acquiring large flow values. However, the QSE cluster continues to shift upward in mass. 
The blockage of the largest flows for the or-captures is a non-equilibrium effect, and the precise 
nuclear mass range where these large flows are localized in depends on the details of the specific 
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expansion profile. For the profile of the bottom row in Figure 10 the largest flows shift in mass 
until the "^"Ca - '^'^Ti - "^^Cr region and produce these isotopes with a mass fraction X ~ 10"^. The 
largest flows shift downward in mass next, while the temperature approaches its minimum value 
at time t = 10"^ s. At this time, the ^*^Si mass fraction reaches its maximum value (bottom right 



panel in Figure 10). The subsequent temperature and density increase until time ? = 1 s relocates 
the largest flows for the or-captures back to the "^^Ca - "^"^Ti - "^^Cr region at the cost of the ^^Si 
abundance. "^"^Ti dominates the final composition with a yield XC^'^Ti) = 0.286. Other isotopes to 
be produced in significant amounts include "^^Ca and ^^Mn. Isotopes near the mass A = 56 are 
underproduced, since their equilibrium abundances are not favored by the QSE formation, and the 
largest flows of a-capture reactions never reach this mass region. 



The non-monotonic profiles of Figure 10 demonstrate cases of final compositions which can- 
not be achieved by monotonic profiles, such as the exponential and power law profiles. Both non- 
monotonic expansions include significant intervals of non-equilibrium nucleosynthesis followed 
by the reformation of equilibrium clusters. This feature results in a final composition dominated 
by silicon group elements which belong to the first family of isotopes, and neutron-rich isotopes of 
the second family near the mass A = 86. Isotopes of the second family near the mass A = 56 are 
underproduced. During monotonic expansions, the QSE cluster size decreases gradually during 
the evolution and the non-equilibrium part of nucleosynthesis is always constrained near freeze- 
out. This feature results in the dominance of isotopes in the second family only, and does not allow 
patterns such as those of Figure [TO} 

The mixing processes during the supernova explosion impact significantly the nucleosynthesis 
mechanisms. The majority of one-dimensional supernova models tend to position the mass-cut 



near the region where the electron fraction begins to decrease below 7^ = 0.5 ( [Woosley et al. 
T973| [Weaver & WoosleyfT993| [Woosley & Weaver) 1 995 1 [Thielemann et al.|l996||Rauscher et al. 
2002[ [Woosley et aL][2002[ [Limongi & Chietti[[2003| [Chiefti & Limongi[[2004[ ). As a result, the 
subsequent supernova shock wave in these models traverses material which has initially symmetric 
or nearly symmetric composition. Our exponential and power-law analysis has demonstrated that 
the radioactivities for A > 60 are not produced in significant amounts during explosions within 
symmetric matter which lack mixing processes. However, there exist types of non-monotonic 



profiles which produce these isotopes even for initially symmetric matter. Figure 1 1 shows two 
such profiles. 



The profile in the top row of Figure 



11 



has extremum values T^'" = 0.157, T^"' = 1.01, 
Pmin = 1.177 X 10^ g cm"^, and p,„ax = 1-13 x 10' g cm"^, and maximizes the ^"Fe yield within our 
data set. During the non-equilibrium part of the evolution the QSE cluster dissolves in multiple 
small-scale clusters. The dominant yields are determined by a-capture reactions within the mass 
range 12 < A < 48, which is a non-equilibrium and profile dependent effect. Yields for masses 
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A > 60 are configured by the interplay of (1) (n, y) reactions out of equilibrium which supply the 
nuclear flow to A > 60 nuclei at the cost of free neutrons, (2) (p, n) reactions out of equilibrium 
which redistribute the nuclear flow among neutron-rich isotopes, and (3) the small-scale chains of 
(p, n) reactions in equilibrium along isobars which collect the majority of nuclear flow available. 
For instance, equilibrium chains along isobars near ^"Fe include ^^Cu to ^^Fe for mass A = 58, 
^''Cu to 5^Co for mass A = 59, '^"Zn to ^°Ni for mass A = 60, ''^Zn to ^^Ni for mass A = 61, and 
^^Zn to ^^Ni for mass A = 62. Although ^''Fe does not participate in the A = 60 equilibrium chain, 
it controls the amount of incoming flow from the Z = 26 isotopic line that is transferred to the 
A = 60 chain through the reactions ^°Fe(p, n)^''Co and ^°Co(p, n)^°Ni. Once the free neutrons are 
depleted, the mass fractions in this mass range are stabilized until freeze-out. The top right panel 
of Figure [IT] shows that ^^Se is also produced. Its mass fraction increase close to freeze-out stems 
from the action of weak reactions (see Section [5]). 

The interplay between (n, y) and (p, n) reactions for initial Ye = 0.5 is another example of 
non-equilibrium nucleosynthesis. The formation of small-scale equilibrium patterns is strongly 
dependent on the particular combinations of temperature and density values during the expansion. 
The profile in the bottom row of Figure [TT] looks similar to the profile in the top row. It has 
extremum values Tg""" = 0.55, T™'-" = 1.41, p„„„ = 3.4 x 10^ g cm-\ and p,„ax = 4.5 x 10^ g cm^^ 
The mass fraction curves in the right column of the Figure have similar trends, but the yields' 
distribution is diff"erent. '^O is underproduced with respect to '^C and ^^Si, and the production 
of ^^Tc and ^^Nb is favored instead of ^°Fe and '^Se. This is an example of the composition's 
dependence on the thermodynamic conditions during the non-equilibrium part of the evolution. It 
is noteworthy though that radioactivities in the mass regime A > 60 are not produced during the 
exponential and power-law profiles for = 0.5, because the large-scale equilibrium patterns for 
most types of freeze-outs (Figure |2]) favor nuclei primarily within the silicon and iron groups up to 
A « 56. 



7. SUMMARY 

We have used parameterized expansion profiles to explore the details of nucleosynthesis trig- 
gered by two mechanisms during freeze-out expansions in core-collapse supemovae. The mass 
layers processed by each of the two mechanisms are separated by a contact discontinuity and do 
not mix. The first mechanism is related to the convection and instabilities within homogeneous 
progenitor matter that is accreted through the supernova shock. The second mechanism is related 
to the impact of the reverse shock on the supersonic wind at the inner regions of the ejecta above the 
proto-neutron star. The exponential and power-law monotonic profiles are nucleosynthesis probes 
for the cases where the supernova shock (or the reverse shock) raises the temperature of the pro- 
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genitor mass layers (or the inner mass layers respectively) above the NSE threshold and subsequent 

mixing processes during the expansion are negligible. Our non-monotonic profiles aim to simulate 
thermodynamic trajectories which are affected by the explosion's asymmetries, instabilities, and 
mixing processes following the passage of the supernova shock through the progenitor mass layers, 
and the effect of the reverse shock on the proto-neutron star material when the temperature does 
not exceed the NSE or large-scale QSE threshold. 

The isotopes produced during the freeze-out expansions are separated in two families. The 
first family of isotopes have a region-divided structure within our parameter space of peak temper- 
atures, peak densities, and initial electron fraction values. Each region in this space is associated 
with a freeze-out type, and its locus depends on the local equilibrium patterns near the isotope 
while the large QSE cluster dissolves. The freeze-out types are characterized by unique equilib- 
rium state transitions (EST) that the QSE cluster, and hence the isotopes of the first family, sustain. 
The mass fraction curves within a region for the isotopes of the first family are similar, and their 
specific profile is shaped by a few critical reactions which differ from isotope to isotope. These 
critical reactions are also responsible for the specific shape of the regions in our parameter space 
and their trends from profile to profile, because they determine whether an EST takes place or not. 
The first family includes the majority of isotopes in the mass range 12 < A < 122. The isotopes 
of the second family are related to nuclei near the magic numbers 28, 50 and 82 and they tend to 
dominate the final composition. These isotopes are produced by maintaining the maximum nuclear 
flows in their locality and they do not sustain any EST for all types of freeze-out. 

The freeze-out types identified within our parameter space involve normal, ff-rich, ap-hch, 
an-nch, the chasm, (p, 7)-leakage and photodisintegration regime. Freeze-out types are classified 
according to the EST that the large-scale QSE cluster sustains, and additional ESTs that isolated 
nuclei sustain in local equilibrium patterns. The local equilibrium patterns are formed once the 
participating nuclei are left outside the QSE cluster and are classified in two general categories. 
The first category results in the mass fraction configuration of nuclei until the iron peak. It involves 
(p, y) reaction chains in equilibrium along isotone lines. If any of the (p, y) reactions along a chain 
breaks equilibrium, then the isotopes of the related isotone line sustain their second EST. Dominant 
yields are localized near the magic number 28. The second category requires significant amounts 
of free neutrons and results in the mass fraction configuration of nuclei beyond the iron peak until 
nuclear masses A ~ 122. It involves (n,y) reaction chains in equilibrium, connected with (p,n) 
reactions in equilibrium. Once the (p, n) reactions break equilibrium the isotopes along the related 
isotopic lines sustain an EST. Dominant yields are localized near the neutron magic numbers 28, 
50, and 82. 

We performed reaction rate sensitivity studies using the exponential and power-law profiles 
and utilized non-monotonic expansion profiles to investigate nucleosynthesis trends of radioactiv- 
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ities from ^^Al to ^^Tc. Once produced, these radioactive isotopes could be transferred to cooler 
mass layers by mixing processes during the explosion, and possibly released to the interstellar 
medium. Contrary to the exponential and power- law profiles, non-monotonic expansions involve 
longer non-equilibrium nucleosynthesis intervals. The production mechanism details are strongly 
dependent on the temperature and density values during the non-equilibrium part of the evolution, 
which implies a dependency on the values of the extremum points for the temperature and density 
trajectories. The non-monotonic expansions demonstrate mass fraction trends and yield distribu- 
tions that cannot be achieved by the exponential and power-law profiles. For instance, there are 
cases where the silicon group yields are larger than the iron group yields, which is not possible for 
monotonic expansion profiles. In addition, the exponential and power-law profiles tend to produce 
^°Fe, ^^Se, ^^Nb, ^^Zr, and ^'Tc only for initially neutron-rich composition, while non-monotonic 
profiles may produce them even for initially symmetric composition. 
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Table 1. 553 
Isotope Nuclear 
Network 



z 


A ■ 


A 


H 


2 


3 


He 


3 


3 


Li 


6 


7 


Be 


7 


9 


B 


8 


II 


C 


II 


14 


N 


12 


15 


O 


14 


19 


F 


17 


21 


Ne 


17 


24 


Na 


19 


27 


Mg 


20 


29 


Al 


22 


31 


Si 


23 


34 


P 


27 


38 


S 


28 


42 


CI 


31 


45 


Ar 


32 


46 


K 


35 


49 


Ca 


36 


49 


So 


40 


51 


Ti 


41 


53 


V 


43 


55 


Cr 


44 


58 


Mn 


46 


61 


Fe 


47 


63 


Co 


50 


65 


Ni 


51 


67 


Cu 


55 


69 


Zn 


57 


72 


Ga 


59 


75 
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Table 1 — Continued 



z 


A ■ 


A 

^max 


Ge 


62 


78 


As 


65 


79 


Se 


67 


83 


Br 


68 


83 


Kr 


69 


87 


Rb 


73 


85 


Sr 


74 


91 


Y 


75 


94 


Zr 


78 


95 


Nb 


82 


97 


Mo 


83 


98 


Tc 


86 


99 


Ru 


89 


99 


Rh 


93 


99 
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Table 2. Isotopes 
of the Second 
Family near the 
Magic Number 28 



z 


A ■ 


A 

^max 


Fe 


56 


57 


Co 


56 


57 


Ni 


56 


62 


Cu 


59 


63 


Zn 


60 


65 


Ga 


63 


67 


Ge 


62 


69 


As 


68 


71 
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Table 3. Nuclear Reactions Relevant to the Synthesis of Radioactivities 



Reaction 


Contribution 


Ye 


Regime 


Profile 




2^Alg(p, 7)2^Si 


Chasm formation, depth, shift 


0.50-0.52 


Chasm 


Both 


2''Alg(a,p)29Si 


Chasm formation, depth, shift 


0.50-0.52 


Chasm 


Both 


26Alg(p,r)"Si 


2nd arc formation/dip 


0.50-0.52 


orp-rich 


Both 


26Alg(p,y)27Si 


3rd arc formation 


0.52 


ap-rich 


Both 


25Al(e-,Ve)'^Mg 


Flow transfer to symmetric nuclei 


0.50-0.52 


orp-rich 


Both 


25Mg(p,y)26Alg 


2nd/3rd arc formation 


0.50-0.52 


orp-rich 


Both 


25Al(p,7)26Si 


2nd/3rd arc formation 


0.50-0.52 


orp-rich 


Both 


24Mg(Q',7)28Si 


Post-arc ascending track 


0.52 


orp-rich 


Power law 


2^P(e-,ye)2^Si 


Flow transfer to symmetric nuclei 


0.50-0.52 


ap-rich 


Exponential 


26Si(e-,Ve)2''AU 


Flow transfer to symmetric nuclei 


0.50-0.52 


ci;;?-rich 


Exponential 


26si(p,r)27p 


2nd arc dip 


0.50-0.52 


a/?-rich 


Exponential 




3rd arc formation 


0.50-0.52 


ap-rich 


Exponential 


2^Si(e-,Ve)2^Al 


Flow transfer to symmetric nuclei 


0.50-0.52 


ap-rich 


Power law 


22Mg(e-,Ve)^^Na 


Flow transfer to symmetric nuclei 


0.50-0.52 


ap-rich 


Power law 


4iCa 


^iCa(p,r)42sc 


Chasm formation, depth, shift 


0.50-0.52 


Chasm 


Both 


^iCa(p,r)42Sc 


3rd arc formation/slope 


0.50-0.52 


or-rich 


Both 


4"Ca(a,r)44Ti 


2nd/3rd arc dip/slope 


0.50-0.52 


Qr;?-rich 


Both 


4iCa(Q;, y)45Ti 


2nd arc formation 


0.48 


tt-rich 


Both 


4«Ca(p,y)4iSc 


3rd arc formation 


0.50-0.52 


ap-rich 


Both 


4iSc(p,r)42Ti 


3rd arc slope 


0.50-0.52 


ap-rich 


Both 


5^Cu(p,7)60Zn 


Flow transfer within QSE cluster 


0.50-0.52 


or-rich 


Both 


3^Ar(Qr,7)4iCa 


2nd arc formation 


0.48 


or-rich 


Both 


^iCa(Qr,p)44Sc 


2nd arc formation 


0.48 


or-rich 


Both 




5^Fe(e",ye)^^Mn 


Flow transfer from symmetric nuclei 


0.50-0.52 


orp-rich 


Both 


52Fe(e-,Ve)^2Mn 


Flow transfer from symmetric nuclei 


0.50-0.52 


orp-rich 


Both 


6«Zn(e-,ye)^''Cu 


Flow transfer from symmetric nuclei 


0.50-0.52 


Qr;?-rich 


Both 


52Mn(p, n)52pe 


Flow transfer from symmetric nuclei 


0.50 


or-rich 


Both 


^2Ti(e-,Ve)^2Sc 


Flow transfer from symmetric nuclei 


0.50-0.52 


ap-rich 


Both 
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Table 3 — Continued 



Reaction 


Contribution 


Ye 


Regime 


Profile 


^^Cr(e-,Ve)^^V 


Flow transfer from symmetric nuclei 


0.50-0.52 


Qfp-rich 


Both 


53Fe(p,7)54Co 


Flow transfer from symmetric nuclei 


0.50-0.52 


ap-vich 


Both 


^3Ni(e-,Ve)"Co 


Flow transfer from symmetric nuclei 


0.50-0.52 


ap-rich 


Both 


^2Fe(p,y)53Co 


Flow transfer from symmetric nuclei 


0.50-0.52 


ap-rich 


Both 




60Fe(p,n)^«Co 


Main depletion reaction 


0.48 


orn-rich 


Both 


^"FeCp, n)^"Co 


2nd arc formation 


0.48 


arn-rich 


Both 




Flow transfer near = 50 


0.48 


a«-rich 


Both 


6"Fe(Q',n)'^3Ni 


2nd arc formation 


0.48 


Q-n-rich 


Both 


^'^Fe(n,7)^"Fe 


2nd arc amplitude 


0.48 


an-rich 


Both 


5^Fe(a,n)6iNi 


2nd arc slope 


0.48 


orn-rich 


Both 


^^Se 


^^As(p,n)'9Se 


2nd arc formation 


0.48 


arn-rich 


Both 


^^Kr(Qr,n)89Sr 


Flow transfer near N = 50 


0.48 


orn-rich 


Both 


^9As(,e- Ve)^^Se 


Flow transfer at freeze-out 


0.48 


an-rich 


Both 


^^Se(n,r)^''Se 


Flow transfer before freeze-out 


0.48 


orn-rich 


Power law 


^^Se(n,7)'^Se 


Flow transfer before freeze-out 


0.48 


ar«-rich 


Power law 




^2Zr(p, nf^m 


2nd arc formation 


0.48 


Q-n-rich 


Both 




Flow transfer near N = 50 


0.48 


an-rich 


Both 




^2Zr(n, 7)9321- 


Flow transfer near N = 50 


0.48 


an-rich 


Power law 




Flow transfer near A' = 50 


0.48 


(wj-rich 


Both 




''^Mo(p,n)''^Tc 


2nd arc formation 


0.48 


an-rich 


Both 


'^'Tc(n,7)'^8Tc 


2nd arc formation 


0.48 


an-rich 


Both 


''^Mo(p,n)''»Tc 


2nd arc formation 


0.48 


Q-n-rich 


Both 


^^Kr(Qr,n)^'^Sr 


Flow transfer near N = 50 


0.48 


an-rich 


Both 
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Table 4. Comparison between the Monotonic and Non-monotonic Expansion Profiles 



Features 



Exponential 



Power-law 



Non-monotonic 



Timescale 

Parameters 
T, p coupled 
QSE evolution 
Dominant yields 
Radioactivities production 



Fixed 
Constant 
1 

No 

Hierarchical 



Dynamic 
Constant 
1 

Yes 
Hierarchical 



Dynamic 
Non-constant 
6 

No 
Periodic 



Second family only Second family only Both families 



< 0.5 



< 0.5 



y, < 0.5 
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Table 5. Final Mass Fractions of Radioactivities for Non-monotonic Expansions 



Ye 


po (g cm 3) 


'rmin 
^9 


^ 9 


Pmin (g cm ^) 


Prruix (g cm ^) 


Mass Fraction 


26^1 


0.48 


10^ 


1.0 


1.2 


2x 10^ 


1 X 10^ 


8 X 10"^ 


0.48 


10^ 


0.04 


1.54 


1.4 X 10" 


1.2 X 10** 


3 X 10-^ 


0.5 


10^ 


0.12 


1 


10^ 


1.5 X 10' 


1.8 X 10^3 


0.5 


10^ 


0.03 


2.32 


2.6 X 10^ 


6.2 X 10^ 


4 X 10-3 


0.5 


10« 


0.44 


0.78 


1.7 X 10^ 


4.4 X 10^ 


10-3 


0.52 


10^ 


0.01 


1.5 


2x 10^ 


10^ 


0.2 


0.52 


10^ 


0.02 


0.9 


2x 10^ 


5x 10' 


4 X 10-3 


^iCa 


0.48 


10^ 


0.18 


2.5 


10^ 


5x 10" 


10-2 


0.48 


10^ 


1 


2 


2x 10^ 


4.3 X 10' 


10-3 


0.5 


10^ 


0.04 


2 


8x 10^ 


4x 10' 


10-3 


0.5 


10^ 


1.7 


2.3 


10^ 


3x 10" 


10-^ 


0.5 


10^ 


1.3 


2.1 


1.4 X 10^ 


6x 10^ 


2 X 10-3 


0.5 


10^ 


0.98 


1.95 


5.7 X 10^ 


2.4 X 10" 


6 X 10-^ 


0.52 


10^ 


1.8 


2.4 


5x 10^ 


5x 10^ 


10-2 


0.52 


10' 


0.08 


0.6 


2x 10^ 


3x 10" 


10-3 


44 


0.48 


10' 


1.4 


1.94 


3.3 X 10^ 


3.4 X 10" 


0.3 


0.48 


10' 


2.3 


1.78 


1.9 X 10^ 


3.8 X 10** 


0.2 


0.5 


10^ 


0.06 


2 


7.1 X 10^ 


8.1 X 10' 


0.4 


0.5 


10^ 


0.99 


1.95 


1.2 X 10^ 


1.1 X 10^ 


0.28 


0.5 


10^ 


1.5 


2.36 


8.8 X 10^ 


3.7 X 10" 


0.18 


0.5 


10^ 


1.77 


1.88 


1.5 X 10^ 


5x 10^ 


0.1 


0.52 


10' 


1.7 


2.3 


10^ 


3x 10" 


0.27 


0.52 


10' 


2.2 


1.86 


3.7 X 10^ 


3x 10^ 


0.19 


53Mn 


0.48 


10' 


1.25 


2 


10^ 


2x 10" 


10-2 


0.48 


10' 


0.2 


3.7 


7.8 X 10" 


10' 


10-2 


0.48 


10' 


0.04 


3.8 


1.5 X 10^ 


3.7 X 10" 


10-2 


0.5 


10^ 


0.02 


2.26 


1.2 X 10^ 


4.3 X 10" 


4.7 X 10-^ 
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Table 5 — Continued 



Ye 


Po (g cm 3) 


f'min 
^9 


'rmax 
^ 9 


Pmin (g cm 3) 


Pmax (g cm 3) 


Mass Fraction 


0.5 


10^ 


2.39 


2.14 


1.2 X 10^ 


1.2 X 10' 


2.7 X 10"^ 


0.5 


10^ 


0.07 


3 


8x 10' 


2.3 X 10^ 


4 X 10-3 


0.5 


10« 


1.67 


2.84 


4.8 X 10^ 


1.2 X 10^ 


1.5 X 10-3 


0.52 


10^ 


2.5 


2.2 


10^ 


2.4 X 10^ 


6 X 10-5 


0.52 


10^ 


0.46 


1 


10^ 


6.8 X 10^ 


5 X 10-5 


60pg 


0.48 


10^ 


0.02 


2.5 


2x 10^ 


3 X 10' 


10-3 


0.48 


10^ 


1.5 


1.3 


10^ 


5.5 X 10^' 


10-3 


0.48 


10^ 


0.04 


0.2 


2x 10^ 


2x 10*^ 


10-5 


0.5 


10^ 


0.04 


0.54 


1.7 X 10^ 


5.8 X 10^ 


4 X 10-^ 


0.5 


10^ 


0.04 


1.14 


5.5 X 10^ 


2.4 X 10^ 


3.7 X 10-6 


0.5 


10« 


0.16 


1 


1.2 X 10^ 


1.1 X 10' 


1.9 X 10-6 




0.48 


10^ 


2 


1.5 


10^ 


2x 10^ 


10-2 


0.48 


10^ 


0.04 


1.5 


4x 10^ 


2x 10^ 


10-3 


0.48 


10^ 


0.04 


0.03 


3 X 10^ 


8x 10^ 


6 X 10-^ 


0.5 


10^ 


0.06 


1.52 


5.6 X 10-^ 


2x 10^ 


10-^ 


0.5 


10^ 


0.06 


0.6 


3.8 X 10^ 


4x 10^ 


6 X 10 5 


0.5 


10^ 


0.34 


0.8 


10^ 


7x 10^ 


5.6 X 10-5 


0.5 


10^ 


0.26 


2 


1.2 X 10^ 


10« 


1.8 X 10-6 




0.48 


10^ 


0.1 


0.5 


7x 10^ 


2x 10' 


2 X 10-2 


0.48 


10^ 


0.1 


1.2 


6.5 X 10^ 


4.2 X 10^ 


2 X 10-2 


0.48 


10^ 


0.08 


2 


10^ 


10' 


10-3 


0.5 


10^ 


0.17 


1.14 


1.1 X 10^ 


9x 10^ 


1.9 X 10-6 


0.5 


10^ 


0.16 


0.38 


2x 10^ 


2x 10' 


1.8 X 10-6 


0.5 


10« 


0.17 


0.29 


1.7 X 10^ 


7.4 X 10^ 


5.8 X 10-5 


0.5 


10^ 


0.55 


1.15 


10^ 


2x 10^ 


2.8 X 10-' 




0.48 


10' 


1.7 


2 


2.5 X 10^ 


10' 


4 X 10-3 


0.48 


10' 


0.07 


1.5 


2.5 X 10^ 


1.5 X 10^ 


3 X 10-3 
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Table 5 — Continued 



Ye 


po (g cm 3) 


'rmin 
^9 


'-pmax 
^ 9 


Pmin (g cm ^) 


Pmax (g cm ^) 


Mass Fraction 


0.5 


10^ 


0.11 


0.4 


5.7 X 10-^ 


2x 10*^ 


7 X 10-' 


0.5 


10^ 


0.11 


1.1 


3.5 X 10^ 


2x 10' 


7 X 10' 


0.5 


10^ 


0.49 


1.11 


1.6 X 10^ 


4x 10^ 


4.7 X 10"' 


9^Tc 


0.48 


10' 


1.9 


2.3 


3.1 X 10^ 


10^ 


5 X 10-^ 


0.48 


10' 


2.7 


2 


3.8 X 10^ 


2x 10' 


4 X 10"^ 


0.48 


10' 


0.1 


2.4 


1.8 X 10^ 


2x 10^ 


2 X 10"^ 


0.5 


10^ 


0.19 


1.7 


1.4 X 10^ 


1.2 X 10' 


7 X 10"^ 


0.5 


10« 


0.19 


0.46 


1.8 X 10^ 


5.8 X 10' 


1.4 X 10^5 


0.5 


10^ 


0.55 


1.41 


3.4 X 10^ 


4.5 X 10^ 


5.7 X 10-' 



-32- 




Time 



Fig. 1. — Schematic temperature or density evolution for the exponential, power-law and non- 
monotonic profiles. Passages through different burning regimes for various peak conditions are 
indicated. The dashed red curve illustrates the impact of variations to the values of the local 
extremum points for the same peak conditions. 
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Peok Temperoture (Tg) 

Fig. 2. — Final mass fractions of the c-chain isotopes (^^Ni omitted) in the peak temperature- 
density plane for the exponential profile at Yg = 0.5. The white colored space corresponds to 
values below the color scale shown. From left to right, the first row corresponds to '^C, '^O and 
^"Ne, the second row corresponds to ^"^Mg, ^^Si and ^^S, the third row corresponds to ^^Ar, "^^Ca 
and '^'^Ti, the fourth row corresponds to '*^Cr, ^^Fe and ^°Zn, and the fifth row corresponds to ^"^Ge, 
^*^Se and '^^Ki. 



-34- 




10 20 30 40 50 

Neutron number (N) 




10 20 50 40 50 

Neutron number (N) 




10 20 30 40 50 

Neutron number (N) 



Fig. 3. — QSE cluster motion in the chart of nuclides for an a-rich freeze-out. The QSE cluster 
remnant condenses near the magic number 28. Each colored line corresponds to a nuclear reaction 
and indicates the level of nuclear flow transferred between the isotopes connected. Normalized 
flows (p are colored black for < < 0.01, navy for 0.01 < (p < 0.05, blue for 0.05 < < 0.1, 
cyan for 0.1 < (p < 0.4, green for 0.4 < < 0.8, red for 0.8 < (p < 1.0, and yellow for 0=1.0. 
Small (p values indicate reactions in equilibrium, while 0=1.0 implies pure one-way nuclear flow 
transfer. 
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Peok Temperature (Tg) 

Fig. 4. — Final mass fractions of isotopes with protons and neutrons near the magic number 28 
(second family) in the peak temperature-density plane for the power-law profile at Ye = 0.52. The 
white colored space corresponds to values below the color scale shown. From left to right, the first 
row corresponds to ^^Ni, ^^Co, ^''Fe and ^°Cu, and the second row corresponds to ^^Zn, ^^Ga, ^^As 
and ^^Ge. 
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Peok Temperoture (Tg) 



Fig. 5. — Final mass fractions of select radioactivities in the peak temperature-density plane for the 
exponential profile at Ye = 0.48. The white colored space corresponds to values below the color 
scale shown. From left to right, the first row corresponds to ^^Al^, '^'Ca, and "^"^Ti, the second row 
corresponds to ^^Mn, ^''Fe, and ^^Se, and the third row corresponds to ^^Nb, ^^Zr, and ^^Tc. The 
a-rich and C3;n-rich freeze-out regions are labeled on the temperature-density plane for '^'Ca. 



-37- 




455789 10 456789 10 



Peok Temperoture (Tg) 

Fig. 6. — Final mass fractions of select radioactivities in the peak temperature-density plane for the 
exponential profile at Ye = 0.52. The white colored space corresponds to values below the color 
scale shown. From left to right, the first row corresponds to ^^Al^, and "^^Ca, and the second row 
corresponds to '^'^Ti and ^''Mn. 
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Fig. 7. — Mass fraction evolution of dominant yields (left panel) and radioactivities (right panel) 
for the transition region between the c-rich and an-rich freeze-outs. Initial conditions are Tg = 9, 
p = 5 X 10^ g cm~^, and Ye = 0.48 for the power-law profile (3304 isotope network). 
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Fig. 8. — QSE cluster motion in the chart of nuchdes for the transition region between or-rich and 
an-rich freeze-outs. The QSE cluster remnant bypasses the isotopes of the second family near 
the magic number 28, and condenses around the corresponding isotopes near the magic number 
50. Each colored line corresponds to a nuclear reaction and indicates the level of nuclear flow 
transferred between the isotopes connected. Normalized flows (p are colored black for < < 
0.01, navy for 0.01 < (p < 0.05, blue for 0.05 < (p < 0.1, cyan for 0.1 < < 0.4, green for 
0.4 < (f> < 0.8, red for 0.8 < (p < 1.0, and yellow for 0=1.0. Small (p values indicate reactions in 
equilibrium, while (p= 1 .0 implies pure one-way nuclear flow transfer. 
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Fig. 9. — Mass fraction schematic for decreasing temperature, where the arc structure is illustrated. 
Such mass fraction profiles are representative of the chasm, a-rich, ap-rich, and orn-rich freeze- 
outs. The first arc (black) describes the mass fraction trend during a large-scale QSE state. The 
second arc (red) describes the mass fraction trend once the isotope is outside the large equilibrium 
cluster, and its trends may be explained by local equilibrium states (local QSE). The blue ascend- 
ing track is denoted as third arc, and is related to a mixture of local QSE and non-equilibrium 
nucleosynthesis. Mass fractions during an a-rich freeze-out have only an ascending track past the 
first arc, while for the chasm region the mass fractions have only one arc. 
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Fig. 10. — Non-monotonic temperature and density trajectories (left column) and key mass fraction 
profiles (right column) for initial conditions = 9, p = 10' g cm"^, and Ye = 0.48. The top row 
is a profile where the ^°Fe yield is approximately maximized, while the profile at the bottom row 
tends to maximize the '^'^Ti yield within our data set. 
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Fig. 1 1 . — Non-monotonic temperature and density trajectories (left column) and key mass fraction 
profiles (right column) for initial conditions = 9, p = 10^ g cm"^, and Ye = 0.50. The top row 
is a profile where the ^°Fe yield is approximately maximized, while the profile at the bottom row 
tends to maximize the ^^Tc yield within our data set. 



